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Abstract 

We study cotunneling in a double junction Coulomb blockade device under 
the influence of time dependent potentials. It is shown that the ac-bias leads to 
photon assisted cotunneling which in some cases may dominate the transport. 
We derive a general non-perturbative expression for the tunneling current in 
the presence of oscillating potentials and give a perturbative expression for 
the photon assisted cotunneling current. 
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I. INTRODUCTION 



Single electron tunneling (SET) devices have attracted much attention during the recent 
years. They can be used to manipulate electrons one by one, and it is the hope that they can 
be applied as e.g. memory cells or as fundamental current standards. However, there are 
still a number of unanswered questions regarding the limitations of single electron devices. 
One such limitation is the so-called cotunneling process where an electron charge tunnels 
through the devices thus giving rise to a leakage current. This process has been studied 
extensively both theoretically^! and experimentally.^ Experiments have confirmed the 
theory of Averin and Nazarovjl who derived the cotunneling current for a double junc- 
tion system at small bias and low temperatures-ll! It has been argued within the present 
understanding of these devices that extremely high accuracy could be achieved leading to 
metrological applications-SBE! However theory and experiments differ by several orders of 
magnitude.!^ At present it not clear what the source of this discrepancy is. The intrinsic 
noise properties of the materials due to charge fluctuations^ is one possible source, but 
the quantum mechanics of the dynamically driven system is also an issue which should be 
considered in order to understand the fundamental limits. 

Therefore an important aspect of single electron devices is their dynamical property 
and their response to an applied high frequency signal. So far theoretical works on single 
electron circuitsSffi have used an adiabatic approximation where it is assumed that the 
static cotunneling rate can be used at every time step. At first sight this approximation seems 
justified when the applied frequency is much lower than A/h, where A is the electrostatic 
energy required to add or remove an electron. As we shall show, this is however not so 
obvious because activated cotunneling is not exponentially suppressed but increases as a 
power law of applied energy. The power law is a consequence of the fact that it is phase- 
space that limits the cotunneling process. 

In this paper we investigate the importance of the photon assisted cotunneling induced 
by time dependent bias and gate voltages in the case of a double junction system in the 
Coulomb blockade regime. Note that this is different from cases previously discussed in the 
literature^0il where the assisting photons are supposed to be supplied from an incoherent 
source, namely the electromagnetic environment connected to a heat bath. The coherent- 
photon-assisted cotunneling is shown to yield an important and in fact for certain parameters 
the dominant contribution to the current. We develop a non-perturbative formalism that 
combines cotunneling and sequential tunneling and thus generalizes previous works which 
dealt with the static case.i'iffil However, our approach is restricted to the weak coupling 
limit and does not include renormalization effects. In the low frequency limit, we separate 
the photon assisted cotunneling in two contributions, one corresponding to an adiabatic 
approximation and one corresponding to higher order processes. In particular, we derive a 
perturbative result for the case of zero dc bias, which we refer to as a cotunneling electron 
pump. 

The paper is organized as follows. In Section [nj we introduce the double junction model. 
The tunnel current is derived in Section |TTT| using scattering theory and several limits are 
considered: weak tunneling, sequential tunneling, the adiabatic limit and the case of small 
voltage and low temperature where a perturbation expansion is applicable. In Section [TV| 
numerical results are given and, finally, conclusions are in Section |V]. Technical details of 
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the time dependent scattering theory are in an appendix. 

II. THE DOUBLE JUNCTION MODEL 

We use the standard tunneling model for transport through a double barrier system with 
Coulomb interactions represented by the electrostatic energies. The model Hamiltonian 
reads 

H = H + H T + H'(t), (1) 

where H = Hi + Hr + Hd + H int and Hl, Hr, and Hp represent the free electron parts in 
the left lead, in the right lead, and in the dot, respectively. The charging energy is included 
in Ho and it is given by 

H mt = E C (N D - n g f = E(N D ), (2) 

where E c is the charging energy, n g is set by the gate voltages and N D is the number of 
electrons in the dot. In Eq. ([I]) is the tunneling Hamiltonian with tunneling matrix 
element Tl(r) connecting the dot to the left (right) lead. The time dependent parts are 
contained in H' and we consider the case where the potentials of the leads and the dot (by 
the gate voltage) are modulated in time 

H'(t)= e a (t)N a} (3) 

a=L,R,D 

where N a is the number operator in section a. 

By a standard unitary transformation (using % = 1): U — exp jz dt' if(f') j , we may 
transform the Hamiltonian into H — ► H + Hx(t), whereby the time dependence is now 
present in the tunneling Hamiltonian 

H T (t)= (W*)dta* + H.c.). (4) 

k,p;/3=L,R 

Here c and d are operators in the leads and in the dot, respectively, and 

u p {t) = exp (i j* dt' [ep(t') - e D (t')]) . (5) 

Without loss of generality we simplify the following discussion by transferring any non-zero 
dc part from the time dependent potential to the static bias or gate voltages. Consequently, 
if e a (t) is periodic so is u(t) and with the same period. 

III. SCATTERING THEORY CALCULATION OF THE TUNNELING CURRENT 

In order to calculate the cotunneling rates we generalize the standard time independent 
scattering theory and obtain an expression for the average transition rate Fj/ between an ini- 
tial and a final state (both eigenstates of the unperturbed Hamiltonian) under the influence 
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of a time dependent perturbation (here Hx{t)). The T- matrix expansion then forms the 
basis for a non-perturbative expression for the current in the limit of small tunneling matrix 
elements. This approach thus generalizes previous scattering theory formulationsaBill! to 
finite frequency and finite temperature. We furthermore utilize a Breit-Wigner-type approx- 
imation which does not include renormalization effects in the strong tunneling regime.^ 
The result for the transition rate is derived in Appendix |A|. After time averaging it reads 
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r fi = 27r^25(E i -E f -u n ) (f\T(n)\i) , (6) 

n 

with u n = 2nn/T and T being the period of the potentials. Here n corresponds to the 
number of photons absorbed. The scattering operator T is 

T{n) = H T (n) + ^ # T (n - m)G (m)H T (m) 

m 

+ 2J Ht{u — to — m')Go(m + to') 

mm 1 

xH T (m)Go(m')H T (m') + ---, (7) 

where the Fourier transforms are defined as F(n) = J dte tuJnt F(t)/TQ. and the unperturbed 
Green's function operator is given by 

G °( m ) = 1 u I ■ • ( 8 ) 

Ei + u m - H + ir] 

The scattering operator T in Eq. @ can be written as 

T(n) = J2 H An - m)G{m, m')H T {m'), (9) 

mm 1 

where we have defined the Green's function operator by the Dyson-like equation 

G(m,mf) = G (m)5 m , m i 

+ J2Go(m)E(m,m")G(m",m'), (10) 

m" 

where 

S(m, to') = H A m ~ m")G (m")H T (m" - m'). (11) 

m" 

Note that we have restricted the summation to even powers of H?, since we shall calculate 
the transition rate of electron transfer through two barries. 



A. Weak tunneling approximation 

When the tunneling is not too strong, the off-diagonal element of the "self-energy" in 
Eq. ([TT|) can be neglected, allowing for a solution of the Green's function. With the notation 
A aa i = (a\A\a'), the Green's function becomes 
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G aa i{m,m) = — — - 12 
[G ,a{ m )\ + s 4 m ) 

where 

Ti a (m) = (a\HT(m — m")\a')Go !a '(m") 

a' ,m" 

x(a'\H T (m" -m')\a). (13) 
Inserting this into Eq. (§), we obtain for the T-matrix 

T if (n) = ^(i\H T (n-m)\a)G a (m)(a\H T (m)\f). (14) 



This formula is equivalent to a Breit-Wigner formula (see e.g. 20), but here it is generalized 
to the case of a time dependent perturbation. 

B. The tunnel current 

Next we calculate the currents through the double junction system, that is the transition 
rate 7 + for transferring an electron from left to right and the rate 7~ for electrons moving 
in the opposite direction. The total current is / = e(7 + — 7~). In the positive direction the 
final state corresponding to an inelastic cotunneling event is given by 

\f + ) = dl c L,pCRyd k ,\i), (15) 

whereas |/~) is obtained by interchanging left and right. The intermediate state denoted by 



a in Eq. flU), can take two values, corresponding to an electron being transfered through the 
left junction or the right junction first followed by an electron tunneling through the other 
junction. We assume that the coherence between the individual tunnelings can be neglected, 
which means that the tunneling uncertainty time %/ max(eV, kT, hui), is much smaller than 
e/I. In this case, the current can be obtained from the transition rate between |/) and \i), 
where |z) is thermally distributed, but with a fixed number of electrons on the island. 

The lifetimes of the intermediate states are given by the imaginary part of S a . The 
real part is neglected since it only gives a small shift of the electrostatic energy. The 



lifetime obtained from Eq. (13) equals the half of the Fermi's Golden Rule transition rate 
for tunneling out of the intermediate state. In the present paper, we restrict ourselves to a 
two state approximation (denoted by N = 0, 1) and therefore the lifetime is given by the 
transition rate from the intermediate state to the initial state. 

By summing over all possible intermediate states while assuming the electronic occupa- 
tion in the initial state to be in thermal equilibrium we finally obtain for tunneling rate in 
the positive direction 



7 + = P(0) J 



00 duj L r°° du) R 



2tt J-oo 2tt 

xT l (uj l + eV L )D LR (u L , uj r )T r (uj r - eV R ), (16) 
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where 

1 — exp(— pu) 

and where P(0) is the probability that the system initially was in the state with N = 0. We 
must include this factor since our derivation for the current assumes that the initial state 
corresponds to the minimum electrostatic energy. 

Above, we have introduced a function D, which describes the propagation of the inter- 
mediate states. It is defined by 

D LR (u 1 ,u 2 ) = 2tcJ2 S ( uj i +^2 -w n ) 

n 

UL{m)u* R (n — m) 



E 



A+ + lji - oj m + iY-UjOi 




+ u* R (m)u L (n - m) . 

A- + u 2 - u m + iT+(uo 2 - w m )/2/ 

Here the first term corresponds to an electron first tunneling through the left junction 
followed by an electron tunneling out of the right junction, while the second term comes 
from the reverse process. In Eq. flUf ) u m is the photon energy in the intermediate state while 
u n is the energy which is absorbed (emitted) during the cotunneling process. A non-zero uj n 
thus represents photon assisted cotunneling. The electrostatic energies in the intermediate 
states are given by 



A ± = E(±l) - £7(0). (19) 



The lifetimes in Eq. dl8f) are 



r±H = r±(u;) 



= En \u a (n)\ 2 T a (uj - uj n =F eV a ] 



(20) 



and is thus the transition rate for an electron to tunnel through one junction in the 
presence of the oscillating fields. These rates are well-known from previous works on photon 
assisted tunneling in Coulomb blockade systemsM@ T a is the tunneling rate without the 
ac field given in Eq. (|T7|). 

The reverse rate is obtained by interchanging left and right. We assume that the bias is 
such that Vl — V/2 and V R = —V/2 (any asymmetry can be absorbed into the voltage on 
the island gate). We then obtain the final expression for the cotunneling current through 
the device as 



e 3 J-oo 2n J-oo 2ix 

/(W! + eV/2)D LR (uj U LU 2 )f(LU 2 + eV/2) 

-f(ui - eV/2)D RL (uj u LU 2 )f(u 2 - eV/2)] . (21) 

Note that a non-zero current can result even for V = if the time dependencies of the 
left and right junctions are different. Below we calculate the "pump" current based on this 
mechanism. 



6 



1. Sequential tunneling limit 



In the limit of small Y the tunneling current reduces to the expression for the sequential 
tunneling described by the familiar master equation approach.il This is seen as follows. 

For small T^, the tunneling rate 7 + only has contributions from the poles of the function 
D. It can be shown that only diagonal terms of D contributions and D thus reduces to a sum 
of delta functions with weights l/r ± (— A T ), and in the two state approximation we keep 
only the terms corresponding to tunneling between N = and N = 1. After integration the 
current then becomes equal to 



-^sequential 

P(0) 



-A+)r+(A+)-r /; (-A+)r+(A- 
r-(-A+) 



(22) 



and from Fermi's Golden Rule it follows that T" 1 
Tij being the rate for tunneling from state N = 
processes. From the master equation we have that P(0) 
sequential tunneling result^! follows. 



(A+) = r i and r - (-A + ) = r 10 , with 
i to N = j including photon assisted 



Tio/ (rio + r i) and the 



C. Harmonic potentials 

In this paper, we focus on harmonically varying fields. We define e a (t) = W a cos(o;ot+0a) 
and find that the functions u become U(L,R){n) = J-n{.0-{L,R)D) and u*{n) = [u(— n)]*, where 
J n are the Bessel functions of the first kind. The arguments am and 4>ld are given by 

a LD = a L + a r>- ^OL L a D cos(0 L - <fi D ), (23a) 
a LD sm(fi LD = a L sm(f) L - a D smcj) D , (23b) 

where ai = Wl/ojq and an = Wd/ujq and with similar relations for L — > R. Using the weak 
tunneling limit, where the energy arguments of the lifetimes in D are approximated by the 
pole values, we obtain after some manipulations 



D LR {u u u 2 ) = 2tt£5(u; 1 + u; 2 -0|£ J m -n{^R D )J m {^LDy m{ - <t ' RD - <l>LD) 

n m 
( 1 1 

L A+ + lo\ - uj m + ir 10 /2 A" - uji + u m + «r 01 /2 



This is the form which is utilized in the numerical results reported in Section [TV]. The D 
function obeys the following relationship 

D rl (uj u u 2} A ± , eV) = D LR {u 2 , u x , A T , -eV). (25) 
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D. Adiabatic limit 



In the low frequency limit where the frequency of the driving fields is much smaller 
than the inverse uncertainty time associated with the cotunneling event, c^o *C A/ft, the 
frequency dependence in the denominator in the expression for D LR , Eq. (|24]) , drops out. If 
it is furthermore assumed that the eV, kT, Wl,r,d "C A, the energy denominators may be 
replaced by A 1 * 1 , and the summation can be performed. This is equivalent to using static 
rekt by Averin and NazaroJ with a time dependent voltage given by eV L (t) — eV R (t) = 
eV + Wl cos(u t + 4>l) — Wr cos(kj t + 4>r)- The current becomes in this approximation 



j(o) _ G L G R h 

adiabatic 127T6^ 



(2*kTf + (eVf + 5(y 2 



X (^ + ^)' (26) 

where W 2 = Wl + W R — 2WlW r cos(<Pl — 4>r)- The two first terms is the result of ref. ||. 
We see that the this formula gives no contribution when no dc bias is applied. 

In the general case, the adiabatic result can be derived from the static cotunneling 
formula, i.e., with Dl R (ui,u)2) setting u — 1, and u = (which is the formula obtained by 
Pasquier et a/.l) but allowing for the gate and lead voltages to vary in time. To next leading 
order in the applied energies eV, kT or Wl,r,d, we obtain after time averaging the formula 

(!) _ G L GrH 2 2 

-'adiabatic _ 247Te 3 RD LD i 



- ,2 



(A+-A-)(A+ + A 
(A+A-) 3 

x[(27rkT) 2 + 3W 2 /A]. (27) 
In the next section, we calculate the lowest order finite frequency correction to Eq. (P?D. 



E. Perturbative result for photon assisted cotunneling current 

Here we derive an perturbative result for the cotunneling current in the case of a "cotun- 
neling pump", i.e. the current without a dc bias. We assume that both the applied frequency 
and the temperature are low enough so that we can expand in powers of the energies. In 
addition, this approximation has the property that one can neglect the broadening effects 
due to the small finite lifetimes, because it is assumed that the energy is so low that the 
real transitions (corresponding to zeroes of the real part of the denominators), do not play 
a role. 

In the case of zero dc bias the cotunneling current is given by Eq. (|2ll) with V = 0. Now 
by using Eq. ( [25] ) and by interchanging oj\ and u 2 in Eq. (|2T| ) the difference D LR — D RL can 
be written as Dl R (A + , A~) — Dl R (A~, A + ). Therefore the pump current vanishes when 
the gate voltage is set to zero. This is expected because in that case there is no preferred 
tunneling direction. 
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Next we expand in powers of the applied frequency uoq whereby the integral in Eq. ([H]) can 
be performed analytically. After summation over frequencies we obtain the final expression 
for the pump current (reinserting H) 

-fpACT.pump = \WrD - W Id\ 

(A+-A-)(A+ + A-) 2 
(A+A-) 3 

x[(2nkT) 2 + 3W 2 /A + (hu ) 2 ]. (28) 

The two first terms is the adiabatic result derived in the previous section. The last term is 
the lowest order quantum correction to the adiabatic limit. 



IV. NUMERICAL RESULTS 



In the following, we compare expression (28) with the current through a Coulomb island 
calculated on the basis of standard master equation approachll, the socalled sequential 
tunneling limit obtained, see Section JIH B 1| . In this calculation we must use the photon 
assisted sequential tunneling rates.0@ In Fig. |l] we show the different current contributions 
as a function of frequency for the case where an ac current is applied to the left lead only. At 
low frequencies the photon assisted cotunneling is much larger than the sequential tunneling 
current. Even with a finite bias voltage (b) this still holds. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have shown that the coherent photon assisted cotunneling is an im- 
portant contribution to the charge transfer processes in single electron devices. It must 
therefore be taken into account when these devices are operated under periodic time de- 
pendent conditions, which is indeed the case in many potential applications. The photon 
assisted contribution that we have pointed out in this paper has not been included previ- 
ously in models for SET systems, therefore it would be interesting to study experimentally 
the photon assisted cotunneling process, e.g., by studying a double junction device under 
the pump conditions, i.e. with zero dc bias but asymmetrically biased. 

However, since the coherent-photon-assisted cotunneling process is important in a small 
parameter space only, it may be difficult to separate the competing tunneling mechanisms 
from the process studied here. In order to observe the photon assisted cotunneling process 
one should thus bias the device as asymmetrically as possible because it is proportional 
to the difference between the ac amplitudes, as can be seen in the perturbative result in 
Eq. (]28D , whereas the usual photon assisted sequential tunneling is additive in the applied 
ac power. Furthermore, since the cotunneling current stems from a higher order tunneling 
processes the resistances should not be too small. In the numerical examples presented 
in Fig. [1], resistances of 0.1 h/e 2 were used. In this case the photon assisted cotunneling 
current dominates for small frequencies but it is important to note that a part of this is 
due to the adiabatic cotunneling current, i.e. it can be explained by a time averaged dc 
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FIG. 1. The photon assisted cotunneling compared to the sequential tunneling current (dashed 
curve) and the perturbative result (dotted curve) derived in Section III E| in the case where only 
the left electrode potential is modulated. The parameters are Wl = 0.5 Ec, Wr = Wd = 0, 
G L = G R = 0.1 e 2 /h, kT = 0.01 E c , and A+ - A - = 0.2E C . Using a typical value for the 
charging energy Ec = -4 meV corresponds to kT = 50 mK and the frequency is scanned up to 40 
GHz in the figure. In (a) there is no applied voltage whereas in (b) eV = 0.1 Ec- 
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cotunneling, giving rise to non-zero time average due to the non-linear bias dependence. 
Since the adiabatic contribution is frequency independent it can in principle be measured in 
the small frequency limit. 

The non-perturbatively expression derived in the this paper generalizes previously de- 
rived cotunneling results for the dc case&HI'El and as in those theories it does not take 
renormalization effects into account. Such effects are important in the strong tunneling 
regime^ and a study of these effects in the time domain is a interesting subject which 
deserves further work. 
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APPENDIX A: TIME DEPENDENT SCATTERING THEORY 

In this appendix, the standard T-matrix approach is generalized to the case of a periodic 
time dependent perturbation, V(t), which is assumed to be turned on adiabatically in the 
distant past, i.e. it includes a factor exp(r]t) where rj = + . We want to calculate the 
probability to find the system at time t in the state |/) given that it started in the initial 
state which is given by the square modulus of the of the overlap, a^j = (f\i(t)). 

First consider the time evolution of an initial state \i) 

\i(t)) = U(t,t )\i(t )), (Al) 

where U(t, t Q ) is the time evolution operator. Expanding in powers of V(t) the nth order 
term becomes (for t — > — oo) 

= (-z)" e -^ot f dh f h ... f tn ~ 1 dt n V I (t 1 )V I (t 2 )---V I (t n )\i), (A2) 

J —oo J —OO J —oo 

where Vj{t) is in the interaction picture. Now define the Fourier transforms V(t) = 
Y^n e ~ tuJ ° nt V{ n ) an d insert these in (A3?) which allows us to perform the integrations and 
obtain 

\i^ n \t)) = e~ iE ' t+rit exp [-iw (mi + • ■ ■ + m n )t] 

{m} 

xG(mi H h m n )V(mi) • ■•G(m n _i + m n )V{m n -i)G Q (m n )V{m n )\i) (A3) 

where G is given by Eq. (g) and u = 2ti/T . 

Taking the square modulus of the overlap aif, we obtain an expression of the form 
(emphasizing only the time dependence) 

l°i/(*)| 2 = £ J2 exp(-iuj [m 1 + --- + m n + m' 1 + --- + m' e ]t)e 27lt --- (A4) 

nl {m}{m'} 
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The time averaged overlap is found by integrating over one period assuming that T <C \jr\. 
(This approximation implies that the frequency of the applied ac signal is much larger 
than the frequency of the tunnel events.) Averaging gives that the sum of the two sets 
of intermediate photon energies must vanish; below we denote the sum of one set as j = 
mi + ■ ■ ■ + me. The tunneling rate is obtained as the time derivative of \a\ 2 and we then 
obtain for the average rate 



2r]e 2vt 



(f\V(j) + E V U ~ m)G (m)V(m) 



( Ei -E f - u j) 2 + V 2 
+ E V(j - m 1 - m 2 )G (m 1 + m 2 )V(m 2 )Go(m 1 )V(mi) ■ ■ 



(A5) 



In the limit i] — > 0, we obtain the result quoted in the main text. 
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